Session  4 


Rotating  Detonation  Engine  Research  at  NRL 


Kailas  Kailasanath 
Naval  Research  Laboratory 


156 


Report  Documentation  Page 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 

VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 

1.  REPORT  DATE 

JUL  2013 

2.  REPORT  TYPE 

N/A 

3.  DATES  COVERED 

4.  TITLE  AND  SUBTITLE 

Rotating  Detonation  Engine  Research  at  NRL 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Naval  Research  Laboratory 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS (ES) 

10.  SPONSOR/MONITOR’S  ACRONYM(S) 

11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release,  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

See  also  ADA593185  2013  International  Workshop  for  Detonations  in  Propulsion  (IWDP  2013)  Held  in 
Tainan,  Taiwan  on  July  26-28, 2013. 

14.  ABSTRACT 

15.  SUBJECT  TERMS 

16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 
ABSTRACT 

SAR 

18.  NUMBER 
OF  PAGES 

14 

19a.  NAME  OF 
RESPONSIBLE  PERSON 

a.  REPORT 

unclassified 

b.  ABSTRACT 

unclassified 

c.  THIS  PAGE 

unclassified 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


ROTATING  DETONATION  ENGINE 
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Focus  is  on  Idealized 
“Engine  or  Device”  Configuration 


Mozzle-End 


Several  groups  around  the  world  are  doing  modeling  and 
experimental  studies  of  RDEs  (Bykovskii,  Wolanski, 
Falempin,  Hayashi,  Schauer,  Yi,  Wang,  Brophy,  Wu,  Clafin, 
Smith,  Tsuboi,  Frolov,  et  al.) 

Recent  RDE  Studies  at  NRL 

Flow-field  description  (2010  International  Combustion 
Symposium) 

Stagnation/back  pressure  effect  (JPC  2010-6880) 

Engine  sizing  effect  (AIAA  201 1-0581) 

Three-dimensionality  effect 
Simulation  of  Specific  rigs 

Various  Fuels  (2012  International  Combustion  Symposium) 
Injection/inflow  effects  (JPC  2011-6044;  ASM  2012-0617,  ASM 
2013-1178) 

Exhaust  flow  (JPC  2012-3943) 

Expansion  Flow  Chemistry 
Preliminary  Fuel-Air  Mixing  studies 

“Propel”  --  More  Efficient  Complex  Configuration  Simulation 
Capability 


How  does  considering  finite 
rate  Chemical  Kinetics  in 
the  expansion  region 
change  RDE  performance. 
Can  this  model  be  used  to 
investigate  pollutant 
formation? 
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BASELINE  SOLUTION 


Baseline  configuration 

-  Stoichiometric  hydrogen-air  RDE  of  Wolanski  and  coworkers* 

-  Geometry 

•  Ratio  of  inlet  nozzle  throat  area  to  wall  is  0.2 

•  14  cm  inner  and  16  cm  outer  diameter  (1cm  thickness) 

•  17.7  cm  axial  length  NOMie-End 

-  Flow  conditions 

•  Hydrogen-air  at  stoichiometric  conditions 

•  Inflow  at  10  atm  stagnation  pressure 

•  300  K  stagnation  temperature 

•  Back  pressure  of  1  atm 


L_. 


1  “  |'  Mlcrornozzlesi 

Head- End 


*Yi,  T-H,  Turangan,  C.,  Lou,  J.,  Wolanski,  P.,  and  Kindracki,  J.  “A Three-Dimensional  Numerical  Study  of  Rotational 
Detonation  in  an  Annular  Chamber,”  AIAA  Paper  2009-634,  Aerospace  Sciences  Meeting,  Orlando,  FL,  2009. 


A-  Detonation  Wave  B-  Oblique  Shock  Wave  D-  Secondary  Shock  Wave 
C-  material  slip  line  between  freshly  detonated  products  and  older  products 

E-  Mixing  region  between  fresh  detonable  mixture  and  detonated  gases 

F-  Region  with  blocked  micro-nozzles,  G-  detonable  mixture  injected  from  micro¬ 
nozzles 
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MODELING  ROTATING  DETONATION  ENGINES 


Geometry  can  be  unrolled  for  2d/3d  simulations 


•  Euler  equations  with  multi-species  mixture  model,  induction  time  parameter  model 

•  Exit  boundary  is  mixture  of  supersonic/subsonic  boundary  with  specified  back  pressure 

•  Typical  small  radial  width  of  combustion  chamber  compared  to  diameter 
allows  us  to  rollout  the  chamber  to  a  two-dimensional  geometry 

•  Domain  extends  to  include  plenum,  injection  plate 

•  Model  does  not  permit  flame  propagation  into  the  mixture  plenum 

•  Inlet  boundary  is  held  at  constant  stagnation  pressure  and  temperature 

Are  these  BC  good  enough  ?  Do  they  capture  potential  back-flow,  pressure  feedback,  swirl  in  exhaust. 

Fundamental  study  to  look  at  Ratio  between  injector  throat  area  and  injection  plate  area  to  see 
how  it  impacts  performance,  feedback 


INJECTION  MODELING 


•  Injector  face 

-  Injector  face  typically  is  a  thick  plate  between  the  combustion 
chamber  and  premixture  or  fuel  and  oxidizer  plenums. 

-  Micro-nozzles  inject  high  pressure  fluid  from  the  plenums  to  the 
combustion  chamber. 

-  The  detonation  wave  typically  travels  adjacent  to  injectors,  so 
injectors  must  be  able  to  withstand  high  heat  flux,  temperatures, 
and  pressures. 

•  Ideal  injector  model 

-  Injector  face  is  specified  as  a  wall  boundary  condition. 

-  Injection  of  the  premixture  is  handled  through  a  source  term  in  the 
density,  momentum,  energy,  and  species  equations. 

-  We  assume  complete  mixing  between  the  injected  mass  and  the 
combustion  chamber  mass  within  one  cell  from  the  injector  face. 

-  Mass,  energy,  and  momentum  of  injected  mixture  is  based  on 
assuming  sonic  conditions  (if  choked)  or  subsonic  conditions  (if 
unchoked)  at  the  exit  plane  of  the  injector. 

-  No  back  flow  is  allowed. 


Relevant  injection 
parameters  are  the  throat 
area  to  the  wall  area 
(A/AJ,  plenum  pressure 
and  temperature. 


-  Valid  for  small  area  ratio  injection  walls. 


Combustion  Chamber 


Premixture  Plenum 
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SLOT  MICRO-INJECTION  SYSTEM 


•  Two-dimensional  simulation 


•  Slot  micro-injector  geometry 

—  Simulation  has  a  total  of  50  micro-injectors,  JE=5.457  mm  spacing 
—  Plate  between  plenum  and  combustion  chamber  is  J=10  mm  thick 
—  Each  injector  has  a  throat  width  of  A=  1.131  mm  for  a=0.2,  2.262  mm  for  a=0.4 


Feedback  into  Mixture  Plenums  in  RDE  s 


50th  Aerospace  Sciences  Meeting,  January  9-12,  2012 


BASELINE  RDE  GEOMETRY 


2 


10 


Injection  , 
Plate  1 

3 


•  8  cm  inner  and  10  cm  outer  diameter 

•  Base  resolution  is  0.02  cm,  1.02  million  cells  for  2D,  51.2  million  for  3D 

•  Injection  plate  has  50  equally  spaced  injectors 

•  Stoichiometric  H2/air  mixture  flowing  into  mixture  plenum 


Exhaust  boundary,  P.  =1  atm 

L L L L L L LJ 


CTTV 

cm 


Buffer 


Combustion  Chamber 


miiiiiiiiiiiiiiiiiiiilllllllli 


Mixture  Plenum 


J  y  J  28.2743  cm  |  |  |  J 

Inflow  boundary,  Pat=4  atm,  T,=300  K 
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INLET  AREA  RATIO  EFFECT 

Pressure  Gain  Combustion 


Set  I. 

Constant  plenum  pressure. 

Area  Ratio 

/plenum  (atm) 

P*. Inlet  (at™) 

/^i.cxii  (atm) 

%  inlet  choked 

0.1 

10 

1.79 

3.92 

87.8% 

0.2 

10 

2.99 

6.66 

90.7% 

0.3 

10 

4.25 

9.39 

94.4% 

0.4 

10 

5.35 

11.91 

94.6% 

0.5 

10 

G.27 

14.08 

94.9% 

0.75 

10 

8.69 

18.38 

95.0% 

Set  2.  Approximately  hold  mass  flow  constant. 


Area  Ratio  F 

plenum  (atm) 

/V. inlet  (atm) 

/Vex it  (atm) 

%  inlet  choked 

0.1 

16.47 

3.13 

6.59 

82.4% 

0.2 

10.00 

2.99 

6.66 

90.7% 

0.3 

6.94 

3.15 

6.72 

72.0% 

0.4 

5.53 

3.21 

6.55 

64.9% 

0.5 

4.71 

3.21 

6.78 

60.3% 

0.75 

3.40 

3.39 

6.21 

26.7% 

*  Inlet  pressure  only  integrated  over  choked  area. 


•  Above  computations  demonstrate  that  for  this  configuration,  to  obtain  pressure  gain  combustion,  we 
must  have  a  large  ratio  of  throat  area  to  wall  area  for  the  injector  face. 

•  As  the  area  ratio  is  increased,  more  feedback  can  be  expected  from  the  combustion  chamber 
detonation  wave  into  the  pre-mixture  plenum  and  compressor. 

•  A  better  understanding  of  how  to  isolate  the  pre -mixture  plenum  from  the  combustion  chamber  is 
essential  in  order  to  achieve  pressure  gain  combustion. 

•  The  simple  mass  addition  injection  model  currently  used  is  insufficient  for  the  higher  area  ratio  cases. 


SLOT  MICRO-INJECTION  SYSTEM 
Animation  to  Illustrate  Pressure  Feedback 


Both  temperature  and 
pressure  show  the  clear 
presence  of  relatively  cool 
jets  at  the  micro-injectors  and 
“dead”  zones  of  higher 
temperature  gas  existing 
between  micro-injector  jets. 

Instantaneous  pressure  field 
also  shows  the  presence  of  a 
trailing  edge  shock  wave  into 
the  mixture  plenum. 

We  have  not  seen  the 
presence  of  black  flow  into 
the  mixture  plenum  in  these 
cases,  although  pressure 
waves  can  be  quite  strong. 
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SLANTED  SLOT  MICRO-INJECTION  SYSTEM 

Geometry 


By  slanting  the  micro-injector  with  respect  to  the  injection  plate  and  main  flow, 
transverse  pressure  waves  propagating  down  the  injector  no  longer  have  a  direct  line 
of  sight  into  the  mixture  plenum. 


•  0  =  0 

Baseline  case  for  comparison. 

•  6  = -20, -40 

Injectors  face  into  the  detonation  wave, 
should  slow  detonation  wave. 

•  0  =  +20, +40 

Provides  maximum  shielding  of  pressure 
pulses  from  the  detonation  wave. 


At  higher  pressures,  we  see  a  reversal  of  the 
detonation  wave  for  high  values  of  0. 


Aw=  5.65  mm 
At  =  2.26  mm  (a=0.4) 
At’  =  2.405  mm  (0=20), 
2.950  mm  (0=40) 


SLANTED  SLOT  MICRO-INJECTION  SYSTEM 
Instantaneous  Pressure  and  Density 


•  Detonation  is  weaker  near  injection  plate 

•  Density  shows  fill  region  inhomogeneity  nicely 
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SLANTED  SLOT  MICRO-INJECTION  SYSTEM 
Instantaneous  Temperature  and  Reactant 


•  Temperature  shows  turbulence  in  detonation,  and  unbumed  region  from  transition. 

•  For  low  pressure  cases,  we  have  broad  transition  region  and  a  small  amount  of 
unbumed  reactants 


SLANTED  SLOT  MICRO-INJECTION  SYSTEM 
Pressure  Trace  in  Mixture  Plenum 


•  5  mm  above  injector  plate  does  not  have  extreme  pressures  seen  20  mm. 

•  Both  overpressures  and  underpressures  exist  in  the  mixture  plenum,  and  are 
as  much  as  +1.4  atm  (+35%)  and  -0.6  atm  (-15%)  from  the  4  atm  stagnation 
pressure  case. 
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SLANTED  SLOT  MICRO-INJECTION  SYSTEM 
Instantaneous  Temperature 


•  +40  degree  angle  no  longer  has  region  of  unbumed  reactants 

•  All  cases  similar  fairly  similar  temperature  field  and  flow- field  features. 


SLANTED  SLOT  MICRO-INJECTION  SYSTEM 
Instantaneous  Density 


•  Both  high  angle  cases  have  very  weak  jets  or  no  jets  at  all  from  injectors. 

•  Still  see  similar  striations  in  the  reactant  density  for  all  cases. 
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SLANTED  SLOT  MICRO-INJECTION  SYSTEM 
Pressure  Trace  in  Mixture  Plenum 


Negative  slant  angles 
give  larger  feedback 
pressures  do  to  how 
disturbances  travel  in 
the  flow-field. 


Positive  slant  angles 
provide  some 
reduction  in  the  over 
and  underpressure 
seen  in  simulations. 
The  overall  values  is 
still  quite  large. 


VARIATIONS  OF  SLOT  MICRO-INJECTION  SYSTEM 

Geometry 


Attempt  to  modify  the  injector  geometry  in  order  to  change  flow  characteristics  in  the 
fill  region  of  the  RDE  and  also  the  amount  of feedback  seen  in  the  mixture  plenum. 


“Cavity  Slot” 


^=5.65  mm 
A  =2.26  mm 
Ac=4.52  mm 
Y  j=7  mm 
Y2=3  mm 


“Nozzle  Slot” 


^=5.65  mm 
A=2.26  mm 
An=5.65  mm 
Yj=4  mm 


“Diode  Slot” 


Aw=5.65  mm 
A  =2. 26  mm 
Ad=4.52  mm 
Yj=7  mm 
Y2=3  mm 


•  Cavity  Slot  -  Operates  by  having  a  small  acoustic  cavity  within  the  injector, 
with  some  energy  of  the  feedback  shock  wave  dissipating  in  cavity 

•  Nozzle  Slot  -  Expansion  near  the  combustion  chamber  to  reduce  size  of  dead 
zones  while  maintaining  same  throat  area. 

•  Diode  Slot  -  Attempts  to  create  a  nozzles  such  that  propagation  of  waves 
upward  is  favorable  compared  to  down  into  mixture  plenum. 
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VARIATIONS  OF  SLOT  MICRO-INJECTION  SYSTEM 
Instantaneous  Temperature  and  Density 


VARIATIONS  OF  SLOT  MICRO-INJECTION  SYSTEM 
Pressure  Feedback  in  Mixture  Plenum 


“Cavity  Slot”  “Nozzle  Slot”  “Diode  Slot” 


•  Cavity  Slot  -  No  substantial  improvement  in  feedback  pressure  compared  to 
baseline  case  and  produces  considerable  turbulence  in  detonation  wave. 

•  Nozzle  Slot  -  Removes  dead  zones  and  jets  in  fill  region  resulting  in 
smoother  fill  zone  and  detonation,  but  is  much  worse  in  terms  of  feedback 
pressure,  since  it  is  more  open  to  the  combustion  chamber. 

•  Diode  Slot  -  No  significant  difference  in  flow-field,  but  slightly  worse  than 
both  the  cavity  slot  and  baseline  case. 


167 


SUMMARY  PERFORMANCE  AND  FEEDBACK 


20  mm  below  plate 

5  mm  below  plate 

Case 

D  (m/s) 

mdot  (kg/s) 

Thrust  (N) 

Isp  (s) 

Pmin(atm) 

Pmax(atm) 

Pmin(atm) 

Pmax(atm) 

Slanted  slot,  angle=0 

1879 

0.530 

594 

4040 

3.49 

5.23 

3.40 

5.36 

Slanted  slot,  angle=-20 

1669 

0.510 

555 

3910 

3.60 

5.45 

3.37 

5.77 

Slanted  slot,  angle=+20 

1913 

0.515 

573 

4010 

3.50 

5.22 

3.37 

5.57 

Slanted  slot,  angle=-40 

1554 

0.492 

532 

3900 

3.58 

5.58 

3.37 

6.05 

Slanted  slot,  angle=+40 

2012 

0.502 

564 

4040 

3.53 

5.24 

3.51 

5.22 

Cavity  slot 

1780 

0.563 

563 

3950 

3.52 

5.26 

3.40 

5.37 

Nozzle  slot 

1810 

0.524 

589 

4050 

3.45 

5.61 

3.44 

6.08 

Diode  slot 

1790 

0.521 

582 

4030 

3.52 

5.46 

3.39 

5.49 

•  Red  and  green  represent  minimum/worst  and  maximum/best  values. 

•  Performance  only  shows  small  variation  (<4%)  between  injectors,  while 
mass  flow  and  thrust  show  larger  variations  (13.7%  and  10.9%) 

•  No  clear  winner  for  reducing  feedback  pressure,  although  nozzle  slot  and 
negative  angle  slanted  slots  are  clearly  the  worst. 


Concluding  Remarks 

•  Rotating  Detonation  Engine  model  developed  at  NRL 

-  Premixed  hydrogen-air  RDE 

-  Low  pressure  ratio  (4),  high  area  ratio  (0.4)  is  regime  of  interest 

-  Baseline  RDE  geometry  with  plenum/injection  plate  simulated 

•  Slanted-slot  micro-injectors 

•  Modifications  of  slot  micro-injectors 
“Cavity”,  “Nozzle”,  and  “Diode”  modifications 

•  Injection  system  results 

-  Flow-field  of  different  slanted-slot  micro-injectors  not  substantially  different 

-  Pressure  feedback  is  minimized  for  positive  angle  slanted-slot  micro-injectors 

-  Alternative  slot  designs  modified  flow-field  in  fill  region  substantially 

-  No  clear  best  design  for  reducing  feedback  pressure  yet,  mixed  bag 

•  Future  work  Further  studies 

•  Evaluate  impact  of  non-premixed  injection 

•  Further  explore  different  isolating  mechanisms  and  geometries 

•  Do  full  “system”  simulation  with  injection  system,  combustion  chamber,  and  exit  nozzle 
design. 
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New  Code  Development  -  “PROPEL” 


Motivation:  Moving  forward,  we  need  to  be  able  to  work  with  a  wider  range  of 
geometries  for  engines  and  experiments,  and  we  need  to  be  able  to  utilize  new 
computational  resources  that  become  available  quickly  while  maintaining  our  ability 
to  simulate  detonation  and  propulsion  concepts  with  the  same  high  fidelity  as  our 
previous/current  work. 

•Based  on  the  Jet-Noise-Reduction  (JENRE)  code  developed  at  NRL. 

•Utilizes  GPU’s,  OpenMP,  Thread-Building-Blocks,  MPI  for  HPC 
•Hybrid  grid  representation 

•  Cartesian 

•  Structured 

•  Unstructured 

•  Combinations 

•Finite  volume  and  finite  element  representations 

•MILES  approach  using  Multi-dimensional  Flux-Corrected-Transport  (FCT) 
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